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Abstract 

Reed-Muller codes are studied for peak power control in multicarrier code-division multiple access 
(MC-CDMA) communication systems. In a coded MC-CDMA system, the information data multiplexed 
from users is encoded by a Reed-Muller subcode and the codeword is fully-loaded to Walsh-Hadamard 
spreading sequences. The polynomial representation of a coded MC-CDMA signal is established for 
theoretical analysis of the peak-to-average power ratio (PAPR). The Reed-Muller subcodes are defined 
in a recursive way by the Boolean functions providing the transmitted MC-CDMA signals with the 
bounded PAPR as well as the error correction capability. A connection between the code rates and the 
maximum PAPR is theoretically investigated in the coded MC-CDMA. Simulation results present the 
statistical evidence that the PAPR of the coded MC-CDMA signal is not only theoretically bounded, but 
also statistically reduced. In particular, the coded MC-CDMA solves the major PAPR problem of uncoded 
MC-CDMA by dramatically reducing its PAPR for the small number of users. Finally, the theoretical 
and statistical studies show that the Reed-Muller subcodes are effective coding schemes for peak power 
control in MC-CDMA with small and moderate numbers of users, subcarriers, and spreading factors. 

Index Terms 

Boolean functions, Multicarrier code-division multiple access (MC-CDMA), Orthogonal frequency- 
division multiplexing (OFDM), Peak-to-average power ratio (PAPR), Reed-Muller codes, Spreading 
sequences, Walsh-Hadamard sequences, Walsh-Hadamard transform. 

I. Introduction 

Multicarrier communications have recently attracted much attention in wireless and mobile applications. 
The orthogonal frequency division multiplexing (OFDM) has been employed as a multiplexing and a 
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multiple access technique in a variety of wireless communication standards such as IEEE802. 1 1 wireless 
LAN Q, IEEE802.16 mobile WiMAX 0, and 3GPP-LTE 0. Also, the multicarrier code-division 
multiple access (MC-CDMA), a combined scheme of OFDM and CDMA 0-[6]], has been proposed 
to enjoy the benefits of OFDM and CDMA by allocating the spread data symbols to subcarriers. The 
popularity of multicarrier communications is mainly due to the robustness to multipath fading channels 
and the efficient hardware implementation employing fast Fourier transform (FFT) techniques. However, 
multicarrier communications have the major drawback of the high peak-to-average power ratio (PAPR) 
of transmitted signals, which may nullify all the potential benefits [7J . 

A number of techniques have been developed for PAPR reduction of OFDM signals. In particular, a 
constructive and theoretical approach is to employ a coding scheme (8] [91] that provides low PAPR and 
good error correction capability for transmitted OFDM signals. The Golay complementary sequences ifTOl , 
which belong to a coset of the first-order Reed-Muller code ifTTl . are a good example of the coding scheme. 
Paterson |[T2l also discussed several coding schemes for PAPR reduction of multicode CDMA. In lfl"3l . 
he summarized the algebraic coding approaches for peak power control in OFDM and multicode CDMA. 
For a summary of the other PAPR reduction techniques for OFDM, see lfl4ll . 

To reduce the PAPR of multicarrier CDMA (MC-CDMA) signals, on the other hand, numerous 
studies have been focused on the power characteristics of spreading sequences. Ochiai and Imai |fT51l 
presented statistical results of the PAPR in downlink MC-CDMA, where multiple users are supported by 
Walsh-Hadamard or Golay complementary spreading sequences. Considering a single user MC-CDMA, 
Popovic fT6l presented the basic criteria for the selection of spreading sequences by studying the crest 
factors (CF) - VPAPR - of various binary and polyphase sequences. Similar studies can be found in 
IfTTl with multiple access interference (MAI) minimization. In MC-CDMA supporting multiple users or 
code channels, the crest factors of various spreading sequences have been compared in |[T8l and |fT9l , 
where the Walsh-Hadamard spreading sequences showed the best PAPR properties, provided that a large 
number of spreading sequences are combined for the transmitted MC-CDMA signals. More studies can 
be found in |[20l — |[23l on the PAPR of various spreading sequences in MC-CDMA. 

If MC-CDMA assigns multiple spreading sequences to a single user, the multicode MC-CDMA can be 
equivalently treated as the spread OFDM (S-OFDM) |[24l — |26l , where a data symbol of the user is spread 
across a set of subcarriers to enjoy frequency diversity. If the number of used spreading sequences is large, 
the Walsh-Hadamard spread OFDM can be viewed as a PAPR reducing scheme ( 181127], compared to a 
conventional OFDM. Also, an error correction code may be applied prior to Walsh-Hadamard spreading 
for improving the error rate performance f26|||2"8l or controlling the peak power |[T8l of S-OFDM. 
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To the best of our knowledge, most of the efforts on PAPR reduction of MC-CDMA and S-OFDM 
have been verified mainly by statistical experiments, not by thorough theoretical analysis. Through the 
experiments, the PAPR of the multicarrier signals has been statistically observed, but it has never been 
addressed whether it is theoretically bounded. In this paper, we propose a binary Reed-Muller coded MC- 
CDMA system and study its PAPR properties. In the coded MC-CDMA, the information data multiplexed 
from users is encoded by a Reed-Muller subcode and the codeword is then fully-loaded to Walsh- 
Hadamard spreading sequences. The coding scheme plays a role of reducing the PAPR of transmitted 
MC-CDMA signals as well as providing the error correction capability. We first establish the polynomial 
representation of a coded MC-CDMA signal for theoretical analysis of the PAPR. A recursive construction 
of Boolean functions is then presented for the Reed-Muller subcodes, where the PAPR of the MC-CDMA 
signal encoded by the subcode is proven to be theoretically bounded. The author of |[29l pointed out that 
the construction is equivalent to Type-Ill sequences in |29l where he made a general and mathematical 
study for Boolean functions with bounded PAPR, not considering the application to MC-CDMA. We also 
discuss a connection between the code rate of the subcode and the maximum PAPR. Simulation results 
show that the PAPR of the coded MC-CDMA signal is not only theoretically bounded, but also statistically 
reduced. In particular, the coded MC-CDMA solves the major PAPR problem of uncoded MC-CDMA by 
dramatically reducing its PAPR for the small number of users. In conclusion, the Reed-Muller codes can 
be effectively utilized for peak power control in MC-CDMA with small and moderate numbers of users, 
subcarriers, and spreading factors. The PAPR properties of the coded MC-CDMA equivalently address 
those of the Reed-Muller coded S-OFDM which supports multiple data from a single user. 

II. System Description 

Throughout this paper, MC-CDMA abbreviates multicarrier CDMA — not multicode CDMA. This 
paper discusses a coded MC-CDMA system employing binary codewords, binary spreading sequences, 
and BPSK modulation. Hence, we focus our description on binary cases. The following notations will 
be used throughout this paper. 

— L is a spreading factor or spreading sequence length. 

— w and W are the actual and the maximum numbers of users supported by a coded MC-CDMA 
system, respectively, where w < W. In the rest of this paper, we will use the context of w users, 
where the w users can be treated as w data bits of a single user in S-OFDM or multicode MC-CDMA. 

— K is the codeword length of a (K, W) code, where K > W. 

— N is the number of information bits that each user transmits in an OFDM symbol. 
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Fig. 1. A general coded MC-CDMA transmitter 



— aW = (ao\«i , • • • ,a>$-i) denotes the iV-bit information of the ith user, < i < w — 1, while 
a n = (a^ , a4 , • • • , a£" l \ 0, • • • ,0) denotes the W-bit uncoded data multiplexed from w users 
and zero-tailed at the nth spreading process, < n < N — 1. Note a"n G {0, 1}. 

— b n = (&n°\ • • • , bn^ ^) denotes the coded output of a n by a (if, code at the nth spreading 
process. Note G {0, 1} for < k < K - 1. 

— d n = (dn \ dn \ ■ ■ ■ , ^) denotes the BPSK modulation output of b„ that experiences the nth 
spreading process. Hence, d^ = (— l) 6 "' G { — 1,+1}. 

— c( fc ) = (c*q\ c[ k \ • • • , c^2 1 ) denotes the L-chip spreading sequence assigned for the A:th coded bit 
of d n , while = (c, , cP"^, • • • ,c\ K l ^) T is a set of the Zth spreading chips across all K spreading 
sequences, where < I < L — l.CisaA'xL orthogonal spreading matrix with L > K, where 
c( fc ) is the kth row vector and c; is the Ith column vector. 

— u n = (un\ Un \ • • • , uw ) denotes the output data of length L of the nth spreading process. 

A. Coded MC-CDMA transmitter 

Figured] illustrates a coded MC-CDMA transmitter proposed in this paper. Assume that w users access 
to the coded MC-CDMA system, where the ith user, < i < w — 1, is actively transmitting the iV-bit 
information aW over an OFDM symbol. The nth information bit of a® is multiplexed across all i's and 



5 



then if w < W, zeros are attached to form the VF-bit uncoded data a n , < n < N — 1. Then, a n is 
encoded by a (K, W) code to generate the iT-bit codeword b n and its BPSK modulation d n . The /cth 
coded bit of d n is then spread by the L-bit spreading sequence c'*', < k < K — 1, where a pair of 
the spreading sequences is mutually orthogonal. At the ?ith spreading process, the spread bits of length 
L are linearly combined over K spreading sequences to produce u n , where each element of u n can take 
an arbitrary value. Obviously, the spreading process is equivalent to a transform of d n by the orthogonal 
spreading matrix C, i.e., u n = d n • C, where c( fc ' is the kth row vector of C. The N blocks of the spread 
data u n of length L experience an N x L block interleaver for frequency diversity, and total N ■ L bits 
are allocated to N ■ L subcarriers by inverse FFT (IFFT). 

The MC-CDMA receiver accomplishes the reverse operation to recover the original information ay' 
for the ith user, where the despreading process is equivalent to a transform by C T , the transpose of C. 

From Figure [T] the baseband transmission signal over an OFDM symbol duration T s is given by 

— N-1L-1K-1 

IEEE 4 fc) c! fc) ^ (Aa+n)t/T * , o < t < t s (i) 

n=0 1=0 k=0 

where j = Note that if the zero-tail processes and the (K, W) encoders are removed from 

Figure [TJ then s{t) with K = w is equivalent to a conventional uncoded MC-CDMA signal in |fT51 . 
The normalization factor a/^T is used in (Q]) to ensure that the average power of s(t) is equal to that of 
the uncoded MC-CDMA signal for w users, which will be shown in Section II-D. 

B. Boolean functions and Reed-Muller codes 

Let x = (xo, • • • ,x m _i) be a binary vector where Xi G {0, 1}, < i < m — 1. A Boolean function 
/(x) is defined by 

2 m -l m-l 

/(x) = f(x , • • • , x m -i) = e ^ IT x% i @) 

i=0 1=0 

where v i € {0, 1} and i\ is obtained by a binary representation of i = Ylh^) 1 *"/2' K k G {0, 1}. Note that 
the addition in a Boolean function is computed modulo-2. In the order of the ith monomial with 
nonzero vi is given by Y^b=o l ^' anc ^ tne highest order of the monomials with nonzero Uj's is called the 
(algebraic) degree of the Boolean function /. 

Associated with a Boolean function /, a binary codeword of length 2 m is defined by 

m—1 

f = (/o, fir - , /2™-i) where fj = /(j'o,ji, ■ ■ ■ ,j m -l), 3 = E J/2 ' 

1=0 
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where ji G {0, 1}. In other words, the associated codeword f of length 2 m is obtained by the Boolean 
function fj while j runs through to 2 m — 1 in the increasing order. 

The rth-order Reed-Muller code R(r, m) is defined by a set of binary codewords of length 2 m where 
each codeword is generated by a Boolean function of degree at most r. In other words, each codeword 
in R(r, m) is the associated codeword f of length 2 m in © where the Boolean function / has the degree 
of at most r. The rth-order Reed-Muller code R(r, m) has the dimension of 1 + (™) + • • • + (™) and 
the minimum Hamming distance of 2 m ~ r . 

For more details on Boolean functions and Reed-Muller codes, see |[30l . 



C. Walsh-Hadamard spreading sequences 

The Walsh-Hadamard matrix is recursively constructed by Hi = [1] and 

H 2 ™-1 H 2 m-1 
H2«"l -Hjm- 



H, 



1 

72 



Then, it is easy to see that 
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where T denotes a transpose and I2™ the 2 m x 2 m identity matrix. §5$ shows that the Walsh-Hadamard 
matrix is symmetric and orthogonal, where the rows (or columns) are orthogonal vectors of length 
2 m , called Walsh-Hadamard sequences. A theoretically defined Walsh-Hadamard matrix lf3TTl has no 



normalization factor in (01). However, we introduce it so that each Walsh-Hadamard sequence has the 
unit energy. 

The Walsh-Hadamard sequences are described by the algebraic structure of Boolean functions and 
the first-order Reed-Muller code. In the Walsh-Hadamard matrix H^™, let h; be the Ith column vector 
of length 2 m , i.e., h, = (h 0;l , h 1>h ■ ■ ■ , h 2m ^ 1>l ) T . Let I = YZ^ 1 ^ and k = YZo* W?, where 
k, ki € {0, 1}. Then, h^ i is given by 



h 



(-1) 



fi(ko,ki,— ,fem-l) 



/nm 



(6) 



where the addition in the exponent is the modulo-2 addition. Without the normalization factor -7=, the 
Ith column vector of length 2 m is a '±1 '-codeword of length 2 m associated with the Boolean function 
fi of m variables. Precisely, /; generates through Q while the row index k runs through to 2 m — 1, 
where is equivalent to a '±l'-codeword in the first-order Reed-Muller code R(l,m). Since each 
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column vector corresponds to a codeword of R(l,m) of the minimum Hamming weight 2 
straightforward that the sum of the column elements is either or 0, i.e., 



m— 1 



i 



k=0 v " fc=0 , 0> 

Example 1: From (01), a 4 x 4 Walsh-Hadamard matrix is given by 

1111 



if I = 
if Z t^O 



it is 



(7) 



H 4 



1 



h. 



1-11-1 
1 1 -1-1 
1-1-1 1 

We have f (k ,fa) = 0, fi(ka,fa) = k , f 2 (k ,fa) = fa, and h(k ,fa) = k + fcj.. Then, it is easily 
checked that each Boolean function /; generates the Ith column vector h;, < I < 3, through (J6J) while 
the row index k = Y2i=o fa^ mns through to 3. Also, (O is true for each column vector. 

The Walsh-Hadamard transform ll30l of a vector g = (go, - " ,92 m -i) is defined by 
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2--1 

= £ fffe -(-l)= 1W % 0</<2 m -l 



fc=0 

where I = Yii^o an d k = Yli=~o~ fa^ 1 - From the algebraic structure of the Walsh-Hadamard matrix 
described above, it is straightforward that the Walsh-Hadamard transform of g is given as 

g = (go,-- - ,?2™-i) =g-H 2 m. 



D. Peak-to- Average Power Ratio (PAPR) 

The peak-to-average power ratio (PAPR) ll32l of a multicarrier signal s(t) is defined by 

FAFR {s(t)) — EUt)l2] 

where T s is an OFDM symbol duration and E[-] denotes the ensemble average. Using the orthogonality 
of spreading sequences, the approach made in ll27ll implies that the average power of the MC-CDMA 
signal s(t) in (Q} is determined by 



iV-l L-l K-l 



2| „(*0|2 



n=0 1=0 k=0 



In particular, if dn £ {— 1 5 +1} and c( fc ) has the unit energy, i.e., Yl 



L-l I Jfc)|2 



E[\s(t)\ 2 ] = — ■ N ■ K = Nw 
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(8) 

1, then ([8]) becomes 
(9) 
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which is equal to the average power of an uncoded MC-CDMA signal where each of w users transmits 
the iV-bit information over an OFDM symbol duration. 

In the following, we define a polynomial S(z) associated with s(t), similar to |[33l . 

Definition 1: In general, the coded MC-CDMA signal s(t) in £Q) has a form of s(t) = EiIo _1 Sie j27rit / T % 
where N ■ L is the number of subcarriers over an OFDM symbol and s, takes an arbitrary value. With 
z = &> 2 ' Kt l T ' , the associated polynomial S(z) is defined by 

NL-l 

S(z) = £ s iZ \ (10) 
i=0 

From © and C[0]), the PAPR of s(t) is translated into 



We establish the polynomial representation of a coded MC-CDMA signal by presenting the associated 
polynomial introduced in Definition Q] For simplicity, we first study the polynomial representation for 
N = 1, where each user transmits a single information bit with a single spreading process in an OFDM 
symbol. The general representation with N > 1 is then discussed. 

A. N = 1 

With N = 1, a coded MC-CDMA signal is denoted by 

S °(*) = \/I ■EEW^- 0<t<T s . (12) 
v i=0 fc=0 

Then, the polynomial representation of so(t) is established by the following theorem. 
Theorem 1: The polynomial Sq(z) associated with so(t) in (fT2l is given by 



where z = (1, z, z 2 , • • • , z L ~ l ) T . In particular, if if = L = 2 m and C is a 2 m x 2 m Walsh-Hadamard 
matrix, i.e., C = H2™, then 



PAPR(s(t)) 



max| 2 | =1 ^(z)) 2 
Nw 



(11) 



III. Polynomial Representation of A Coded MC-CDMA Signal 




(13) 




(14) 



where d is the Walsh-Hadamard transform of d . 



Proof. In CEU>, let s (t) = Zt=o ufe^ 11 ^ where 4° = ^/f • Y,k=o cf ] . Then, 



|.d„.c. 



u = (%\--- ,Uq l) ) = 

With z = ei 2nt / Ts , the associated polynomial is then given by So(z) = Yld=Q u ^o zl 
derives ([T3T ). If C = H 2 m, then (fT4l ) is immediate. 



(15) 

uo • z, which 

□ 



Corollary 1: With N = 1, if a coded MC-CDMA signal s (i) has the PAPR of at most P, then 

|5o(^)| 2 < wP 

from (fTTb . Also, we have from ( TT3T ) 

|d -C-z| 2 = — • |5 (z)| 2 < KP. 



where z = (1, z, z , • • • , z ) with \z\ = 1. 

In (fT5T ). if C = H2>«, the spread output uo is the Walsh-Hadamard transform of do with a scaling 
factor in a Walsh-Hadamard spread MC-CDMA, where K = L = 2 m . A similar property has been 
noticed in a spreading process of multicode CDMA lfT2ll . In addition, H 2 >™ z has the following structure. 



Lemma 1: Let H 2 >« be a 2 m x 2 m Walsh-Hadamard matrix and z = (1, z, z 2 



V 2 m -1\T 



) J . Then, 
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Thus, (fT6l ) is true for 



where zi = 
• z 

where z = 



m = 1. Assume (1161) also holds for m = A; — 1, i.e., 

Go (00) " " ' ) 0fc-2, 00, • • • ) 0fc-2) 
Gl(0o,- -- , 0fc-2)^O> • • • 

Il2fc-i " Z^ 

G 2 fe-i_i(0O) ■ ■ ■ , 0fe-2i ^0, ■ ■ ■ 
(1, z, z 2 , • • • , z 2t 1_1 ) T . From the recursive construction of H 2 t, we have 



(18) 



1 
71 
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— H2*-i 

,2 fc -l\T 
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ofc-i 
2 • Zi 
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• H^k-l " Zx 

• H2*-i • z\ 



0fc-l • H 2 fc-i 
0fc-i " H 2 fc-i 



zi 



(19) 



(l,z,z 2 ,--- ,z 2fc " 1 ) T = (zf,z 2fc " -zf) T . Thus, from (HQ and COO), 

0fc-l -Go(0o,--- ,0fc-2,0O, 
0fc-l • Gi(0o, • • • , 0fc-2) 00, 



H^fc • z 



t7fc-2j 

) 0fc-2) 



0fc-l • ^fe-i-l^O) • • • ) 0fc-2, ^0) ' ' ' ) 0fc-2) 

0fc-l • Go(0O) • • • ) 4>k-2,0o, • • • ) ^fe-2) 
0fc-l • Ci(0o, • • • , 0fc-2, 00, • • • , 0fc-2) 



• G r 2fc-i_l(0O, 



^fc-2, 00, • • • , 0fc-2) 



Gq(4>0, ■ ■ ■ , 4>k-li 00, • • • , 0fc-l) 

G ! i(0o,--- , 0fc-i, 0o, • • • , 0fc-i) 

^2^-1 (00, ' ' ' , 0fc-l, 00, - - - , 0fc-l) 



From (fTTT ) and (l20l ). (fT6l ) is true by induction. 



(20) 

□ 



Remark 1: In Definition 6 of ll34l , Parker and Tellambura defined G m , a set of normalized complex 
sequences of length 2 m by a tensor product, where |0 t | 2 + \9t\ 2 = 2. In fact, H2»> • z in Lemma Q] is a 



special case of G m with 



1+z 2 ' 



and t = ±=#^, < i < m - 1. 

V2 



B. N > 1 



In (fT2l . replacing by d^f* leads us to s n (t) and its associated polynomial S n (z), i.e., 

s n(*) = V¥'EE d n )c P )ei27rZt/Ts ' S n (z) = <M-d n -C* 1 0<n<iV-l (21) 



i=0 fc=0 

»L-X\T 



where z = (l,z,z , ••• ,2: ) . Obviously, s n (i) is also a coded MC-CDMA signal where d n is 
loaded with a single spreading process over an OFDM symbol. In particular, if C = H 2 m, then the 



1 1 



(23) 



spreading process is equivalent to the Walsh-Hadamard transform (WHT), which enables the efficient 
implementation of spreading and despreading processes. 

Theorem 2: With S n (z) in (1211 ). the associated polynomial of a coded MC-CDMA signal s(t) in £[]) 
is determined by 

N-1 

S(z) = £ S n (z N ) ■ z n . (22) 

n=0 

In other words, S{z) is a polynomial obtained by interleaving 5 n (z)'s for < n < N — 1. 
Proo/ In ©, 

— 7V-1L-1K-1 AT-1 
n=0 1=0 k=0 n=0 

where s n {t) = M ■ d^ c\ k) e^ m+n ^ . 

» i=0 fc=0 

For a given n, s n (i) is a signal assigned to the (Nl + n)th subcarriers while I runs through to L — 1. 
From d23l) . it is straightforward that the associated polynomial of J n (t) is given by 

S n (z) = ^-d n -C-z^.z n 

where = (1, z 27V , • • • , z( L " 1 ) Ar ) T . Compared to (EB, = 5„(z Ar ) • z n , and thus the 

associated polynomial S(z) is 

JV-l JV-l 

S(z) = S n (z) = S n(* N ) ■ z n - 

n=0 n=0 

□ 

Remark 2: Since s n (t) is a coded MC-CDMA signal with a single spreading process over an OFDM 
symbol, Corollary Q] is also valid for S n (z). Precisely, if s n (t) has the PAPR of at most P, then |5* n (z)| 2 < 
wP and |d n • C • z| 2 < KP for < n < N - 1, where z = (1, z, z 2 , • • • , z L ' x ) T with \z\ = 1. 

Using its associated polynomial S(z) in Theorem |2j we determine the PAPR bound of a coded MC- 
CDMA signal s(t) with N > 1. 

Theorem 3: In (1211 ). assume the maximum PAPR of s n (t) is P, i.e., maxo< n <jv-i PAPR(s n (i)) = P. 
Then, the coded MC-CDMA signal s(t) in © has the PAPR of at most NP, i.e., 

PAPR(s(i)) < NP. 
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Proof. From the associated polynomial S(z) in dZZl . 

|S(z)| 2 = \S (z N ) + Si(^) ■ z + ■ ■ ■ + SW-i(^) ■ z 



iV-ll 2 



N-l N-l 



n=0 n=0 




where \z\ = 1. From Remark|2l maxo< n <AT-i PAPR(s n (i)) = P implies IS^r^)! < wP for every n. 
Thus, |5"(z)| 2 < N 2 ■ wP. Therefore, the PAPR of s(t) is bounded by 



Although the proof is straightforward and the bound seems not so tight, Theorem [3] gives us an insight 
that the maximum PAPR of coded MC-CDMA signals increases as each user transmits more data bits 
(N) in an OFDM symbol. Therefore, N should be as small as possible to remove the probability that 
the MC-CDMA signal has the high PAPR. 



In this section, we develop a variety of subcodes of R(r, m) for a (K, W) coding scheme of a coded 
MC-CDMA in Figure [T] where the codeword of length K = 2 m is associated with a Boolean function 
of degree r. We assume that the K-bit codeword is fully-loaded to all the available Walsh-Hadamard 
spreading sequences of length 2 m , so K = L = 2 m . We analyze the PAPR properties of the fully-loaded, 
Reed-Muller coded, and Walsh-Hadamard spread MC-CDMA signals. First of all, we study the PAPR 
for N = 1. Then, the PAPR for N > 1 is investigated. 

In the MC-CDMA system with N = 1, d = ((—l) b< " > , (— l) 6 ^' , • • • ) (_l)6? m_1) ) where b = 
^o^'"'' ' ^) i s a codeword of a (2 m ,W) Reed-Muller subcode. In this section, we denote 
b = (b , h, • ■ • , &2 m -i) for simplicity. 

A. The first-order Reed-Muller code 

Let bo = (&0)&i>"" >^2 m -i) be a codeword of the first-order Reed-Muller code = R(l,m). 

When it is employed as a (K, W) coding scheme in a coded MC-CDMA, the dimension is W = m + 1 
and the codeword length is K = 2 m . Each codeword is associated with a Boolean function of 



PAPR(s(t)) 



max| z | =1 \S(z)\ 2 N 2 ■ wP 



= NP. 



Nw ~ Nw 



□ 



IV. Reed-Muller Codes for MC-CDMA 



h(x , ■ ■ ■ 



) = E V i X i + 



e, Vi,ee{0, 1} 



(24) 
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where the addition is computed modulo-2. The PAPR of the MC-CDMA signal encoded by a codeword 
in is determined in the following. 

Theorem 4: With K = L = 2 m and N = 1, let s (t) be a Walsh-Hadamard spread MC-CDMA signal 
encoded by b = (6 , h, • • • , & 2 ™-i) G = 12(1, m). Then, the PAPR of s (t) is 

PAPR(s (t)) = 1. 

Proof. From Theorem [TJ the associated polynomial of so(t) is given by So(z) = ■J^£ • do • z where 
d is the Walsh-Hadamard transform of d , and z = (1, z, z 2 , • • • ,z 2m ~ 1 ) T . By definition, d = 

(do , di , • • • ,d 2 ™_i) where 



2 m — 1 m— 1 m— 1 

£ ( _ 1) fc+EE; 1 W. where I = £ J,*, k = £ fc,2*. 



k=0 i=0 i=0 




In (El), 6 fc = 6i(fc , • • • , fcm-i), < A; < 2 m - 1. From ©, 

2--1 2--1 

. /nm Z-^i * ' /om ^ ' 

v fc=0 

if i = « 
otherwise 

where « = XJIo 1 for § iven ^' s - Therefore, 5 (z) = • d • z = ■ • z v = ±y/w ■ z v . 

For any v, the PAPR of so(t) is therefore 

PAPR( So( t)) = maX '^ l|g ° (z) ' 2 = ^1 = 1. 

□ 

From Theorem [4] we see that the first-order Reed-Muller code is a simple and effective coding scheme 
that provides the uniform power for the coded MC-CDMA signals. However, it has a relatively low code 
rate l?i = which vanishes as the code length increases. Therefore, we need to develop high-rate 

coding schemes at the expense of the PAPR increases. 



B. Recursive construction 

From a seed pair of codes, we present how to recursively construct a new code using the associated 
Boolean functions. We also analyze the PAPR of the coded MC-CDMA signals. 

Theorem 5: Let / and g be Boolean functions of (m — 1) variables, where by G F and h g E Q are 
the codewords of length 2 m ~ l associated with / and g, respectively. Assume that the code rate of each 
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code is Rf and R g , respectively. Let Sf(t) and s g (t) be the coded MC-CDMA signals encoded by by 
and b g , respectively, each of which has a form of ([T2]) where K = L = 2 m_1 and N = 1. Assume that 
each signal has the PAPR of at most P. 

Consider a Boolean function b of m variables defined by 



6(x , • • • ,X m _i) = (1 + %-l) • /(x , • • • ,X m _ 2 ) + £m-l • 5(^0, • • • ,Xm-2)- 



(25) 



Then, a codeword bo of length 2 m associated with b has the code rate R = —hi — -. Let so(t) be a coded 
MC-CDMA signal of CU) encoded by b , where K = L = 2 m and N = 1. Then, the PAPR of s (t) is 

PAPR(s (t)) < IP. 



Proof. Obviously, | F\ = 2 2 



and \Q\ 



R », respectively. Thus, the number of codewords bo 



is \P\ -\g\ = 2 2m l i R f+ R 9) and the code rate is R 



■(Rf+Rg) 



££+^.. in particular, if F = Q, 



then we keep the code rate R = Rf = R g while the codeword length doubles. 

Let df and d g be the BPSK modulation outputs of length 2 m ~ 1 from hf and b g , respectively. From 
d25l ), it is straightforward that bo = (by | b 5 ) and do = (df \ d g ), where '|' denotes a concatenation. 
Then, the associated polynomial Sq(z) is determined by 



S (z) 



■w 



— • do • H9" • z 

2m u 



W 

(df I d g ) 





2m I d g ) 



1 


U 2 r 


n-l 


H 2m 


-1 




Zl 


71 


H 2 , 


n-l 


H 2 „ 






2 m-l 
Z • Zi 


H 2 


m-1 • 


Zi • 


^m- 1 






H 2 


m-1 • 


zi • 


Om-1 







(26) 



(l,2;,z 2 , 
J=-. Let Bf(z) -- 



where z 

_ ™ _ 1 

1-z 



-1\T 



)% Zi = (l,Z,Z 2 ,--- ,Z 2 -f, Vm _i 



-l\T 



1+Z ^ , and 6» m _i = 
Tpr • dj • H 2 n«-i • zi and B g (z) = y^Ir • d 9 • H 2 ™-i • zi. Then, (l26l ) becomes 

S (z) = Bf(z) ■ (j) m -l + B g (z) ■ Om-l- 



Replacing the Boolean function g by g + 1 leads us to the change of the above associated polynomial to 
S' (z), i.e., 

S' (z) = Bf(z) ■ (j) m -i - B g (z) ■ 6 m -i. 



Then, 



\So(z)\ 2 + \S> Q (z)\ 2 = 2.{\B f (z)\'< 



-i\ 2 + \B g (z. 



1 2 • |#m-l| 2 ) 



(27) 
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If s f (t) and s g (t) have the PAPR of at most P, then Corollary Q] implies |d/ • H 2 ™-i • zi| 2 < 2 m " 1 P 
and \dg ■ i • zi| < 2 m ~ P, respectively. Thus, 

\B f (z)\ 2 <^-2 m - 1 -P=^-P, \B g (z)\ 2 < ^ ■ P. (28) 

Therefore, from ((27J and (l28l 

|S (z)| 2 + |S W| 2 < 2 • max (^/(z)! 2 , \B g {z)\ 2 ) ■ {\<j) m -i\ 2 + |# m -i| 2 ) 

= 2- — P-2 = w2P 
2 

where |</> m _i| 2 + |# m _i| 2 = 2 from the definition of 4> m -i and Thus, the PAPR of so(t) is 

PAPR(s (i)) = 111 V 71 < 2P. 

w 

□ 

The recursive construction of a Boolean function has been originally discussed in |[T2l for the PAPR 
of multicode CDMA. In Theorem \5\ we showed that the construction of (I25T ) also provides the bounded 
PAPR for multicarrier CDMA. In general, if there is a seed code B^ 1 ^ of length 2 m_1 and size \B^™ X ^|, 
then we can construct a new code B^ of length 2 m and size |£>r"^| = l-S^j" 1 ^ 2 by concatenating a 
pair of codewords from the seed. If the PAPR of each coded MC-CDMA signal for the seed B^. r \ ^ is 
at most P, then each coded MC-CDMA signal encoded by the new code Br provides the PAPR of at 
most 2P. If each codeword in B r _ x is associated with a Boolean function of degree at most r — 1, then 
B r m ^ is a subcode of R(r, m) defined by a Boolean function of degree r, where the minimum Hamming 
distance of B r is at least 2 m ~ T '. 

Construction Q] summarizes a recursive code construction for the application to MC-CDMA. 

Construction 1: For positive integers r and m, 2 < r < m, let b\(xo, ■ ■ ■ ,x m - r ) = Y^bLft ViXi + e 
and bi(xo, • • • , x m - r ) = Y^I!=q v^Xi + e', where Vi, v' i: e, e' G {0, 1}. Starting with b\ and b[, the Boolean 
function b r (xQ, ■ ■ ■ , a; m _i) of degree r is constructed by the (r — 1) successive recursions of 

(29) 

while f runs through 2 to r. In (1291 ), the Boolean function 6j ! _ 1 has the same form as b?-i, but may have 
different coefficients. Let bo = (bo,bi,--- ,62 m -i) be a codeword of B r C R(r, m), associated with 
a Boolean function of b r . Then, Br has total 2( m ~ r+2 )' 2r 1 codewords through the (r — 1) recursions. 
In a Walsh-Hadamard spread MC-CDMA with K = L = 2 m and A r = 1, the PAPR of s (i) encoded 
by bo is at most 2 r_1 from Theorems |4] and [5] 
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Finally, the code parameters of B^ are summarized as follows. 



Dimension W = 2 r • (m — r + 2) and code length K = 2 m , 

('" 
— 



Code rate i? r = 2 "~ 1 ' ( ™- r+2) = 2 r " 1 - m • (m - r + 2), 



— Minimum Hamming distance > 2 m r , 

- PAPR(s (t)) < 2 r ~ 1 . 

First of all, we present a specific code example B^ of length 2 m through a single recursion, where 
a pair of codewords in R(l, m — 1) is employed as the seed. 

Construction 2: Let bo = (&o 5 , • • • , &2 m -i) be a codeword of B^ C R(2,m) that is associated 
with a Boolean function b 2 defined by 

1 m— 2 

&2O0, • • • ,x m ~i) = ^ ViXi + x m _i • ^2 v 'i x i + e > «i,«i,e G {0, 1}. (30) 

In a Walsh-Hadamard spread MC-CDMA with K = L = 2 m and TV = 1, the code parameters including 
the PAPR of a coded MC-CDMA signal so(t) encoded by bo are summarized as follows. 

— Dimension W = 2m and code length K = 2 m , 

- Code rate R 2 = ^r, 

— Minimum Hamming distance > 2 m ~ 2 , 

- PAPR(soO)) < 2. 

With a single recursion (r = r = 2) in Construction [T] (l30b is straightforward by 

(m-2 \ /m-2 \ 

^ VtXt + e + x m _i • ^ ufxj + e' 

m— 2 m— 2 



^ ^iX, + (e + e') • x m _i + x m -i • ^ (uj + i> •') • + e 



i=0 i=0 
m— 1 m— 2 

= ^ + x m _i • ^ f-Xj + e 

i=0 i=0 

where u m -i = e + e' and u' = u j + 



Remark 3: By generalizing d30l to 

m— 1 7 
i=0 i=0 

we obtain a code QB™ associated with g 2 . Obviously, the coding scheme B^ in Construction [2] is a 
special case of QB { 2 m) for 7 = m — 1. It is not so hard to prove that a coded MC-CDMA signal so(t) 
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encoded by a codeword associated with g 2 has the PAPR of at most 2 for any 7. While 7 runs through 

Item- 1, we have 2 m+1 • (2 m ~ 1 + 2 m ~ 2 - 1 + 2 m ~ 3 - 1 + • • • + 2 - 1) = 2 m+1 • (2 m - m) distinct 

codewords in QB^ 1 ' , more than the number of codewords in B^ l \ If we compare the code rates of 

QB 2 an d , however, the code rate difference of 

log 2 (2 m+1 • (2 m - m)) - 2m _ log 2 (2 m - m) + 1 - m 1 
2"j 2 m ^ 2 m 

is very small and approaches to as m increases. Meanwhile, the encoding and the decoding complexities 

of the generalized coding scheme QB2 are obviously larger than those of £> 2 "^ . For the little contribution 

to the code rate and the increase of the complexities from QB^\ we therefore consider B^ as our 

coding scheme for low PAPR. 

By employing the coding scheme B^, the Walsh-Hadamard spread MC-CDMA system with K = 
L = 2 m and N = 1 is able to support maximum 2m users or 2m information bits from a single user in 
an OFDM symbol, providing the PAPR of at most 2. 

Construction 3: Let a Boolean function 63 be defined by 

m— 1 m— 3 m— 2 m— 3 

b 3 (x , ■■■ , X m -l) = ^ V i X i + ^m-l^m-2 • ^ v[xi + X m -l ■ ^ v'-Xi + X m - 2 • ^ ufiCi + e (31) 
i=0 i=0 i=0 i=0 

where «j, v^v", v"' , e € {0, 1}. Let bo = (bo, b\, • • • , &2 m -i) be a codeword of ,63"^ C R(3, m) that is 
associated with 63. In a Walsh-Hadamard spread MC-CDMA with K = L = 2 m and N = 1, the code 
parameters including the PAPR of a coded MC-CDMA signal so(t) encoded by bo are summarized as 

— Dimension W = 4(m — 1) and code length K = 2 m , 

- Code rate i? 3 = f^, 

— Minimum Hamming distance > 2 m ~ 3 , 

- PAPR(s (t)) < 4. 

Similar to Construction |2l the Boolean function 63 is immediate from the twice recursions of ( f25T ) for 
r = 2 and 3 with r = 3. With the coding scheme , the Walsh-Hadamard spread MC-CDMA system 
with K = L = 2 m and N = 1 is able to support maximum (4m — 4) users or (4m — 4) information bits 
from a single user in an OFDM symbol, providing the PAPR of at most 4. 

Table U lists the parameters of B^ of length 2 m from Construction [TJ where 1 < r < 5 and m > r. It 
elucidates a connection between the code rates of B^ and the maximum PAPR of MC-CDMA signals 
encoded by B^\ Obviously, we obtain a high rate Reed-Muller subcode at the cost of high PAPR for 
the coded MC-CDMA signals. 
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TABLE I 

The parameters of Bl m) of length 2 m for some r's 



Code 


Dimension W 


Degree r 


Minimum 


Maximum PAPR 






of a Boolean function 


Hamming distance 


(for K = L = 2 m and N = 1) 


Si 


m + 1 


1 


> 2 m_1 


1 


B 2 


2m 


2 


> 2 m ~ 2 


2 


B 3 


4(m - 1) 


3 


> 2 m ~ 3 


4 


B 4 


8(m - 2) 


4 


> 2 m ~ 4 


8 


B 5 


16(m - 3) 


5 


> 2 m " 5 


16 



Remark 4: In Construction [TJ if r = m, then R r = 1 and PAPR(s (*)) < 2" 1 " 1 in ^ m) . In other 
words, if we apply a trivial mth-order Reed-Muller code of length 2 m and code rate 1.0, or equivalently 
apply a Reed-Muller mapping to the information bits of length 2 m , then the PAPR of the corresponding 
MC-CDMA signals is bounded by 2 m ~ 1 . The Reed-Muller mapping may be useful for the applications of 
small spreading factors, e.g., m = 2, 3, 4, requiring low PAPR and no rate loss. However, if the spreading 
factor is large, we may have the problems of the high PAPR and the large complexity of demapping at 
the receiver. 

C. Golay complementary sequences 

Golay complementary sequences iflOl provide the bounded PAPR < 2 for transmitted OFDM signals 
when they are employed as a coding scheme in an OFDM system. From ifTTll . it is well known that 
each binary Golay complementary sequence of length 2 m is equivalent to the second-order coset of the 
first-order Reed-Muller code, where the associated Boolean function is defined by 

m— 2 m— 1 

b c (xo, ■ ■ ■ ) = XI X n(i) x n(i+1) + ^2 ViXi + e (32) 

i=0 i=0 

where it is a permutation in {0, 1, • • • , m — 1} and V{, e G {0, 1}. We now apply a binary code for 
coded MC-CDMA signals, where each codeword is a Golay complementary sequence defined by b c . 

Construction 4: Let bo = (feo^i) - "" )^2 m -i) be a codeword of C R(2,m) that is associated 
with a Boolean function b c in (l32l . Then, its code parameters are summarized as ifTTi 

— Dimension W = m + log 2 (m!) and code length = 2 m , 

- Code rate R c = m+1 f/ ml) , 



Minimum Hamming distance > 2 



m-2 
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TABLE II 

The maximum PAPR of Walsh-Hadamard spread MC-CDMA signals encoded by Golay complementary 

SEQUENCES OF LENGTH 2 m (K = L = 2 m , N = 1) 



m 


Code length 


Theoretical maximum (Theorem [6} 


Actual maximum (Numerical experiments) 


3 


8 


4 


1.9654 (< 2) 


4 


16 


4 


1.9998 (< 2) 


5 


32 


8 


3.2184 (< 4) 


6 


64 


8 


3.8826 (< 4) 


7 


128 


16 


3.9930 (< 4) 


8 


256 


16 


5.8964 (< 8) 



Through a linear unitary transform (LUT), Parker and Tellambura |f34l have implicitly investigated the 
PAPR of MC-CDMA signals encoded by the Golay complementary sequences. Therefore, Theorem [6] is 
immediate from the application of their work in a Walsh-Hadamard spread MC-CDMA. 

Theorem 6: O In a Walsh-Hadamard spread MC-CDMA with K = L = 2 m and N = 1, let s (t) 
be the coded MC-CDMA signal encoded by a codeword bo G Be m Construction 01 Then, 

PAPR(soW) < 2 m ~ L H (33) 

Proof. From Lemma [T] T = H2*™ • z is a special case of a linear unitary transform (LUT) described in 
Theorem 6 of j34l . where (f>t = and Ot = ■ Theorem Q] implies that so(t) is equivalent to a 

unitary transform of a modulated Golay complementary sequence do through T. Therefore, the bounded 
PAPR of (O is obvious from Corollary 6 in O. □ 

Remark 5: Numerical experiments revealed that the actual maximum PAPR of MC-CDMA signals 
encoded by Golay complementary sequences is smaller than the upper bound predicted by Theorem [6] 
This may be the case because the PAPR bound of Theorem [6] has been established by general <p t and 6 t 
with the requirement of \4> t \ 2 + \9 t \ 2 = 2 (33. However, the Walsh-Hadamard spread MC-CDMA has the 

113* 1—7 2 ' 

special values of cp-t = and 6 t = , respectively, which may require the tighter bound on the 

PAPR than (l33l . Table [TT] shows the PAPR comparison between the theoretical bound and the numerical 
results. From the numerical results, we conjecture PAPR(so(i)) < 2^" 5~J, where the proof is left open. 

Table [In] compares the code rates of B^™\ B^ m \ and B^ for several m's. We see from the table 
that B^ is a proper coding scheme for a Walsh-Hadamard spread MC-CDMA with a small spreading 
factor, providing the bounded PAPR < 4 and the acceptable code rates. 
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TABLE III 

The code rates of various codes of length 2 



m 


Code length 


B 2 


Ba 


B c 






(Construction |2) 


(Construction [5} 


(Construction f4} 


3 


8 


0.75 


1.0 


0.6981 


4 


16 


0.5 


0.75 


0.5366 


5 


32 


0.3125 


0.5 


0.3721 


6 


64 


0.1875 


0.3125 


0.2421 


7 


128 


0.1094 


0.1875 


0.1508 



D. Encoding and decoding 

In Figured] the Reed-Muller subcode introduced in Construction Q] is applied at each spreading process 
for encoding the information from a single or multiple users in the Reed-Muller coded and Walsh- 
Hadamard spread MC-CDMA transmitter (K = L = 2 m ). Precisely, a (2 m ,W) Reed-Muller subcode 
B r m ^ encodes a VF-bit input data a n = (an , • • • , On l \ 0, • • • ,0) at the nth spreading process, < 
n < N — 1, to produce a codeword b n of length 2 m , which goes through Walsh-Hadamard spreading, 
interleaving, and IFFT in the sequel. 

In the encoding process, the codeword b n is obtained by 

b n = a n • GH 

where is the W x 2 m generator matrix of Br m \ where W = 2 r ~ 1 (m — r + 2). The recursion of 



Boolean functions in d29| ) equivalently derives the recursion of generator matrices of 



G 



(m—r+f) 



, (m— r+f— 1) 



'f-1 







G 





(m— r+f— 1) 



2 < f < r 



where = (0, • • • , 0) of length 2 m ~ r+f - 1 and G 
elementary row operations, it is equivalent to 



f-i 

(rn-r+f) ■ f-l 



is a 2 r " 1 (m - r + 2) x 2 



G 



(m—r+r) 



G(m—r+f—l) 
f-1 



G 
G 



(m—r+f—l) 
f-1 

(m—r+f—l) 
f-1 



2 < r < r. 



(34) 

~ r+T ' matrix. By 
(35) 



While r runs through 2 to r, the generator matrix G r "^ is constructed by the (r — 1) recursions of (1341 or 



(I35I ). where the initial matrix G^ 



(m— r+1) 



is the (m — r + 2) x 2 



m— r+1 



generator matrix ofR(l,m — r + 1) 
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given by 





1111 


1111 • 


• 1111 


1111 




1 






ni m 

Ulul 


UXU-L 


UXU-L 




x 


r+l) _ 


0011 


0011 • 


• 0011 


0011 




Xl 




0000 


1111 • 


• 0000 


1111 




X2 




0000 


0000 • 


• 1111 


1111 




Xm— r 



For the notations Xj's of the generator matrix of the first-order Reed-Muller codes, see 1301 . 

In particular, we are able to determine and G3 directly from the Boolean expressions in (l30l) 

and (f3TT >. respectively. The generator matrices of £>2 and B3 are 2m x 2 m and (4m - 4) x 2 m matrices, 
respectively. With each of length 2 m , we have 

1 

x 



G 



(m) 



1 

XQ 

X m -1 
X m -l x 

Xm— lX m _ 2 



,(m) 



(36) 



Xm — 1 

X m -l x 

Xm— lX m _2 
X m -2 X 

Xm— 2X m — 3 
Xm— 1X774—2X0 

^m— lX m — 2X m _3 

Note that and G 3 "^ in (l36l ) have the different orders of rows with those generated by the recursions 



of (El or ([35]). In this paper, we use G^ m) and G 3 mJ in @fj. 

Example 2: Let m = 3. Then, the generator matrices of £>2 and B% are 6x8 and 8x8 matrices, 



,(m) 
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respectively. 



.(3) 



1111 


1111 




1 


0101 


0101 






0011 


0011 




Xl 


0000 


1111 




x 2 


0000 


0101 




x 2 x 


0000 


0011 




x 2 xi 



G 



(3) 



1111 


1111 




1 


0101 


0101 




xo 


0011 


0011 




Xl 


0000 


1111 




x 2 


0000 


0101 




x 2 x 


0000 


0011 




X 2 Xi 


0001 


0001 




XlX 


0000 


0001 




X 2 XiX 



We briefly introduce some decoding techniques for Reed-Muller subcodes. The first-order Reed-Muller 
code can be decoded by the Fast Hadamard Transform (FHT) technique described in ll35l . In general, 
the rth-order Reed-Muller code is decoded by the Reed decoding algorithm ||36l . In particular, if we 
consider 3% or 63 as a supercode of the union of cosets of R(l,m), then we can accomplish the 
soft decision decoding by removing each possible coset representative from the received codeword and 
then applying the FHT (37) CCD- 

For the encoding and decoding of Golay complementary sequences, see ifTTI . 



E. PAPR of coded MC-CDMA signals with N > 1 

In what follows, we discuss the PAPR of coded MC-CDMA signals in a general case of N > 1. We 
restrict our attention to a Walsh-Hadamard spread MC-CDMA employing B^ or B^ that provides 
the acceptable code rate as well as the low PAPR for the coded MC-CDMA signals. We show that the 
maximum PAPR depends on the actual number of users supported by the MC-CDMA. 

Theorem 7: Assume that B^ is employed in a Walsh-Hadamard spread MC-CDMA system in 
Figure [l] where K = L = 2 rn . The maximum PAPR of the coded MC-CDMA signal s(t) is then 
determined by 

N, if 1 < w < m + 1, 
2N, if m + 2<w < 2m. 

Similarly, if the system employs B^" -1 , then the maximum PAPR of s(t) is 

N, if 1 < w < m + 1, 
PAPR(s(t)) < < 2N, ifm + 2<w<2m, 

AN, if 2m + 1 < w < 4m - 4. 




(37) 



(38) 
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Proof. In Theorem [3l it is easy to see that P = maxo< n <Ar-i PAPR(s n (t)) = maxd PAPR(so(i)), 
where so(t) and s n (t) are given in ([12]) and (|2Tb . respectively. Therefore, P = 1,2, and 4 when sj"^, 
^ m) , and are employed as the coding scheme, respectively. In of (l36l) . if u> < m + 1, the 

first (to + 1) rows participate in the encoding process, while the other rows are ignored by zero tailing. 
Since a linear combination of the first (m + 1) rows generates a codeword of b[ , it is obvious that if 
w < m + 1, then PAPR(s(i)) < N from Theorems [3] and 01 If m + 2 < w < 2m, on the other hand, 
PAPR(s(t)) < 2N from Construction [2] and Theorem [3] Therefore, (I37T ) is true for Z?2 • Similar to this 
approach, d38T ) is also true for 23g from the generator matrix G 3 "^ of (|36l ) and Theorem [3] □ 

In general, if is employed as the coding scheme, the Walsh-Hadamard spread MC-CDMA signals 
have the PAPR of at most N • 2 r ~ 1 from Construction [T] and Theorem [3] However, Theorem |7] is not true 
for the MC-CDMA signals if the generator matrix Gr"^ is recursively constructed by (1341 or (|35l ). We 
need to reorder the rows of to achieve the maximum PAPR depending on the number of actual 

users as in Theorem [7] 

Remark 6: In Section IV-B, we developed various Reed-Muller subcodes of length 2 m to control the 
peak power of MC-CDMA signals with 2 m subcarriers in a systematic way. In fact, we may employ the 
coding scheme for a codeword of length 2 m N, where N = 2 h , which encodes a = (ao | ai | • • • | ajv_i), 
a concatenation of N uncoded data block. Then, the codeword covers the entire 2 m N subcarriers to control 
the peak power and to ultimately reduce the maximum PAPR of the coded MC-CDMA signal. In this 
case, however, the code rate may be dramatically reduced for such a long codeword because the coding 
scheme is a subcode of the Reed-Muller code. This also enlightens a connection between the code rates 
and the maximum PAPR of coded MC-CDMA signals. 

Remark 7: Virtually treating a single user's data as multiple users' one, a coded MC-CDMA system 
can be considered as an equivalent spread OFDM E51 , where the single user's data is spread across a set 
of subcarriers to enjoy frequency diversity. By applying the coding schemes introduced in this section, 
the spread OFDM additionally has the benefits of low PAPR and good error correction capability. 

V. Simulation Results and Discussions 

This section provides simulation results to confirm our theoretical analysis and presents some discus- 
sions on statistical results of PAPR of MC-CDMA signals. The PAPR properties of a Reed-Muller coded 
and Walsh-Hadamard spread MC-CDMA system are compared to those of a pair of uncoded systems. In 
the uncoded systems, the one employs Walsh-Hadamard (WH) spreading sequences, while the other uses 
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Golay complementary (GC) spreading sequences each of which forms a row of a recursively constructed 
Golay complementary spreading matrix |fi~5l . In our coded MC-CDMA, we employ 63 as the coding 
scheme which we believe is a good coding solution providing the acceptable code rates, the moderate 
complexity, and the low PAPR for the coded MC-CDMA signals. 

For a fair comparison, we assume that all the MC-CDMA systems transmit the same number of infor- 
mation bits in an OFDM symbol from w active users. If the uncoded systems transmit Nw information 
bits in an OFDM symbol, our coded system of code rate R then needs to transmit coded bits for 
the transmission of Nw information bits. Therefore, while the uncoded ones have NL subcarriers, the 
coded system needs to use subcarriers in an OFDM symbol, where L is a spreading factor used in 
the uncoded systems. In the following, our simulations employ -83 of code rate R = 0.5 or 1, where 
the coded MC-CDMA uses 2NL or NL subcarriers in an OFDM symbol by employing the spreading 
sequences of length 2L or L. 

In our simulations, we measure the discrete-time PAPR ll38l of each MC-CDMA signal from the IDFT 
(Inverse discrete Fourier transform) of the oversampling factor 8. Also, we statistically measure the PAPR 
over N s = 5 x 10 5 OFDM symbols for NwN s randomly generated information bits. 

A. Code rate R = 0.5 

In Figures [2] and [3) w users access to each MC-CDMA system to transmit N = 4 information bits 
per each user in an OFDM symbol. The uncoded MC-CDMA systems use the spreading factor L = 16 
and NL = 64 subcarriers, where maximum 16 users are supported, i.e., 1 < w < 16. As the coded 

(5) 

MC-CDMA system also needs to support up to 16 users, we choose a (32, 16) code B- A as its coding 
scheme, where W = 16, K = 32, and R = 0.5. Thus, our coded MC-CDMA uses the spreading factor 
2L = 32 and 2A^L = 128 subcarriers to transmit Nw information bits in an OFDM symbol from the w 
active users. In the coded MC-CDMA system, each 32-bit codeword is fully-loaded to all the available 
32 Walsh-Hadamard spreading sequences regardless of w. On the other hand, the uncoded MC-CDMA 
systems assign the w spreading sequences to w users on demand, so they are fully-loaded only if w = 16. 

Figure [2] shows the complementary cumulative distribution functions (CCDF) of Pr(PAPR > Ao) 
of each MC-CDMA signal for w = 8 and 16. It reveals that the coded MC-CDMA is superior to the 
others when the number of active users is small. Precisely, if w = 8, it reduces the PAPR Ao achieving 
Pr(PAPR > Ao) = 10~ 3 by more than 2 dB, compared to the uncoded systems. Moreover, Theorem |7] 
ensures that there exists no coded MC-CDMA signal with PAPR > 9 dB for w = 8, which implies 
that the coded MC-CDMA also outperforms the uncoded ones in theoretical aspects. If w = 16, the 
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Fig. 2. PAPR performance of MC-CDMA systems. All the MC-CDMA systems transmit Nw = 32 or 64 information bits in 
an OFDM symbol. The code rate of the coded MC-CDMA is 0.5. 



coded MC-CDMA has almost the same PAPR A as the uncoded Walsh-Hadamard spread MC-CDMA 
for achieving Pr(PAPR > Ao) = 10~ 3 . Even in this case, it is theoretically guaranteed that no coded 
MC-CDMA signal has PAPR > 12 dB, which may not be true in the uncoded systems. Figure [2] also 
shows that most of the coded MC-CDMA signals in the statistical experiments have much smaller PAPR 
than the theoretical maximum predicted by Theorem [7] 

Figure |3] displays the PAPR A of each MC-CDMA achieving Pr(PAPR > A ) = 10~ 3 according to 
the number of active users w, 1 < w < 16. It is well known |[T5l that the uncoded Walsh-Hadamard 
(WH) spread MC-CDMA shows the high PAPR when the number of active users is small. The PAPR then 
decreases as the number of users increases. On the other hand, the uncoded Golay complementary (GC) 
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Fig. 3. Relationship between PAPR Ao and active number of users for MC-CDMA systems where Pr(PAPR > Ao) = 10 
The code rate of the coded MC-CDMA is 0.5. 



spread MC-CDMA has the low PAPR for the small number of users. However, the PAPR gets higher 
than that of the uncoded Walsh-Hadamard spread MC-CDMA as the number of users increases. Figure [3] 
shows that the coded MC-CDMA is a good alternative to the two uncoded systems by providing the 
smallest PAPR Ao for almost all user numbers. Moreover, Theorem [7] assures that the maximum PAPR 
of the coded system is theoretically limited to 6 dB for 1 < w < 6, 9 dB for 7 < w < 10, and 12 dB for 
11 < w < 16, respectively, where Figure [3]provides the numerical evidences. Therefore, it is theoretically 
guaranteed that there exists no coded MC-CDMA signal with the PAPR higher than the maximum values 
for each user, which is not generally true in the uncoded systems. The theoretical and statistical results 
show that the coded MC-CDMA dramatically reduces its PAPR for the small number of users, which 
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effectively solves the high PAPR problem in the uncoded MC-CDMA. Figure [3] also shows that if the 
number of active users is large (w > 11), the statistical PAPR Ao is much smaller than the theoretical 
maximum predicted by Theorem [7] As a result, we claim that the coded MC-CDMA provides the best 
statistical and theoretical solution for PAPR reduction for any number of users. 

B. Code rate R = 1.0 

In Figures |4] and [5J the MC-CDMA systems support w users where each user transmits = 8 
information bits in an OFDM symbol. The uncoded MC-CDMA systems use the spreading factor L = 8 
and NL = 64 subcarriers to transmit Nw information bits in an OFDM symbol. To support up to 8 users, 
the coded MC-CDMA system employs B^\ a (8,8) code with K = W = 8 and R = 1. In this case, 
B^ is used as a mapping scheme mentioned in Remark 01 In our MC-CDMA, each 8-bit uncoded data 
a n , < n < 7, is transformed by the Reed-Muller mapping scheme for PAPR reduction. Thus, it uses 
the spreading factor 8 and 64 subcarriers to transmit Nw data bits in an OFDM symbol, which is the 
same as the uncoded systems. However, note that each 8-bit codeword of our MC-CDMA is fully-loaded 
to all the available spreading sequences of length 8 regardless of w. 

Figure @] shows the results of Pr(PAPR > A ) of each MC-CDMA for w = 4 and 8. We observed 
that if w = 4, our Reed-Muller mapped MC-CDMA system reduces the PAPR A by about 3 dB to 
achieve Pr(PAPR > A ) = 10~ 3 , compared to the uncoded Walsh-Hadamard (WH) spread MC-CDMA 
employing the same spreading factor and the same number of subcarriers. Moreover, Theorem [7J ensures 
that our system has no signal with PAPR > 9 dB for w = 4. Note that Theorem |7J determines the 
maximum PAPR of 9 dB (if w < 4), 12 dB (if 5 < w < 6), and 15 dB (if 7 < w < 8), respectively. 
Thus, even if the statistical PAPR property of our MC-CDMA is almost identical to that of the uncoded 
Walsh-Hadamard spread MC-CDMA for w = 8, our system has no probability of signals with the PAPR 
higher than 15 dB, which is however unclear in the uncoded systems. 

Similar to Figure [3 Figure [5] displays the PAPR A of each MC-CDMA achieving Pr(PAPR > A ) = 
10~ 3 according to the number of users w, 1 < w < 8. It shows that our MC-CDMA provides the smallest 
PAPR Ao for almost all user numbers. Also, it numerically confirms that the theoretical maximums of 
PAPR in Theorem |7J hold for each user number. Figure [5] shows that if the number of active users is 
large, the statistical PAPR Ao is much smaller than the theoretical maximum predicted by Theorem [7J 
For the small number of users, on the other hand, the coded MC-CDMA solves the high PAPR problem 
of the uncoded MC-CDMA by dramatically reducing its PAPR. Along with Figure [3j the coded MC- 
CDMA can be the best statistical and theoretical solution for PAPR reduction for any number of active 
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Fig. 4. PAPR performance of MC-CDMA systems. All the MC-CDMA systems transmit Nw = 32 or 64 information bits in 
an OFDM symbol employing the spreading factor 8 and 64 subcarriers. The code rate of the coded MC-CDMA is 1.0. 

users regardless of its code rate. A drawback of the Reed-Muller mapped MC-CDMA is that it could be 
employed only for a small spreading factor due to the demapping complexity at the receiver. 

VI. Conclusion 

This paper has presented a coded MC-CDMA system where the information data is encoded by 
a Reed-Muller subcode for the sake of PAPR reduction. In the system, the codeword is then fully- 
loaded to Walsh-Hadamard spreading sequences, where the spreading and the despreading processes are 
efficiently implemented by the Walsh-Hadamard transform (WHT). We have established the polynomial 
representation of a coded MC-CDMA signal for theoretical analysis of the PAPR. We have then developed 
a recursive construction of the Reed-Muller subcodes which provide the transmitted MC-CDMA signals 
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Fig. 5. Relationship between PAPR Ao and active number of users for MC-CDMA systems where Pr(PAPR > Ao) = 10 
All the MC-CDMA systems use the spreading factor 8 and 64 subcarriers. The code rate of the coded MC-CDMA is 1.0. 



with the bounded PAPR as well as the error correction capability. We have also investigated a theoretical 
connection between the code rates and the maximum PAPR in the coded MC-CDMA. Simulation results 
showed that the PAPR of the coded MC-CDMA signal is not only theoretically bounded, but also 
statistically reduced by the Reed-Muller coding schemes. In particular, it turned out that the coded 
MC-CDMA could solve the PAPR problem of uncoded MC-CDMA by dramatically reducing its PAPR 
for the small number of users. Finally, the theoretical and statistical studies exhibited that the Reed-Muller 
subcodes are effective coding schemes for peak power control in MC-CDMA with small and moderate 
numbers of users, subcarriers, and spreading factors. We believe this work gives us theoretical insights 
for PAPR reduction of MC-CDMA and S-OFDM by means of an error correction coding. 



30 



References 

[I] IEEE Standard 802.1 1-2007, IEEE Standard for Information Technology - Local and Metropolitan Area Networks - Specific 
Requirements, Part 1 1 - Wireless LAN Medium Access Control (MAC) and Physical Layer (PHY) Specifications. 

[2] IEEE802.16e-2005, IEEE Standard for Local and Metropolitan Area Networks, Part 16 - Air Interface for Fixed and Mobile 

Broadband Wireless Access Systems, Amendment 2: Physical and Medium Access Control Layers for Combined Fixed and 

Mobile Operation in Licensed Bands. 
[3] 3GPP TS 36.211, v. 8.3.0, Technical Specification Group Radio Access Network; Evolved Universal Terrestrial Radio Access 

(E-UTRA); Physical Channels and Modulation (Release 8). 
[4] N. Yee, J. P. Linnartz, and G. Fettweis, "Multi-carrier CDMA in indoor wireless radio networks," in Proc. of IEEE PIMRC, 

pp. 109-113, Sep. 1993. 

[5] K. Fazel, S. Kaiser, and M. Schnell, "A flexible and high performance celluar mobile communications systems based on 

orthogonal multicarrer SSMA," Wireless Personal Communications, vol. 2, pp. 121-144, 1995. 
[6] S. Hara and R. Prasad, "Overview of multicarrier-CDMA," IEEE Commun. Mag., vol. 35, no. 12, pp. 126-133, Dec. 1997. 
[7] K. G. Paterson, "Generalized Reed-Muller codes and power control in OFDM modulation," IEEE Trans. Inform. Theory, 

vol. 46, no. 1, pp. 104-120, Jan. 2000. 
[8] A. E. Jones, T. A. Wilkinson, and S. K. Barton, "Block coding scheme for reduction of peak to mean envelope power ratio 

of multicarrier communication schemes," Electron. Lett., vol. 30, pp. 2098-2099, 1994. 
[9] R. D. J. van Nee, "OFDM codes for peak-to-average power reduction and error correction," in Proc. of IEEE GLOBECOM, 

London, U.K., pp. 740-744, 1996. 
[10] M. J. E. Golay, "Complementary series," IRE Trans. Inform. Theory, vol. IT-7, pp. 82-87, 1961. 

[II] J. A. Davis and J. Jedwab, "Peak-to-mean power control for OFDM, Golay complementary sequences, and Reed-Muller 
codes," IEEE Trans. Inform. Theory, vol. 45, no. 7, pp. 2397-2417, Nov. 1999. 

[12] K. G. Paterson, "On codes with low peak-to-average power ratio for multicode CDMA," IEEE Trans. Inform. Theory, 

vol. 50, no. 3, pp. 550-559, Mar. 2004. 
[13] K. G. Paterson, "Sequences for OFDM and multi-code CDMA: Two problems in algebraic coding theory," HPL-2001-146, 

Hewlett-Packard Laboratories, 2001. 
[14] S. H. Han and J. H. Lee, "An overview of peak-to-average power ratio reduction techniques for multicarrier transmission," 

IEEE Wireless Communications, pp. 56-65, 2005. 
[15] H. Ochiai and H. Imai, "OFDM-CDMA with peak power reduction based on the spreading sequences," IEEE Int. Conf. 

on Commun. (ICC), vol. 3, pp. 1299-1303, 1998. 
[16] B. M. Popovic, "Spreading sequences for multicarrier CDMA systems," IEEE Trans. Commun., vol. 47, no. 6, pp. 918-926, 

1999. 

[17] S. Nobilet, J. F. Hlard, and D. Mottier, "Spreading sequences for uplink and downlink MC-CDMA: PAPR and MAI 

minimaization," European Trans. Telecommun., vol. 13, pp. 465-473, Sept./Oct. 2002. 
[18] B. -J. Choi, E.-L. Kuan, and L. Hanzo, "Crest-factor study of MC-CDMA and OFDM", IEEE VTC'99(Fall), pp. 233-237, 

1999. 

[19] B. -J. Choi and L. Hanzo, "Crest factors of complementary-sequence-based multicode MC-CDMA sginals," IEEE Trans. 

on Wireless Commun., vol. 2, no. 6, pp. 1114-1119, Nov. 2003. 
[20] E. Pogossova, K. Egiazarian, and J. Astola, "Spreading sequences for downlink MC-CDMA transmission," IEEE 

VTC2004(Fall), vol. 7, pp. 4859-4863, 2004. 



31 



[21] L. Yang and E. Alsusa, "Dynamic code-allocation based PAPR reduction technique for MC-CDMA systems," in Proc. of 

WCNC 2007, pp. 628-633, 2007. 
[22] Q. Shi and Q. T. Zhang, "Symmetry-embedded spreading sequences for multicarrier CDMA," IEEE Trans, on Wireless 

Commun., vol. 6, no. 10, pp. 3534-3539, Oct. 2007. 
[23] D. A. Wiegandt, Z. Wu, and C. R. Nassar, "High-throughput, high-performance OFDM via pseudo-orthogonal carrier 

interferometry spreading codes," IEEE Trans. Commun., vol. 51, no. 7, pp. 1123-1134, July 2003. 
[24] S. Kaiser, "On the performance of different detection techniques for OFDM-CDMA in fading channels," in Proc. of IEEE 

GLOBECOM, pp. 2059-2063, 1995. 
[25] D. Merouane, P. Loubaton, and M. de Courville, "Spread OFDM performance with MMSE equalization," in Proc. of IEEE 

Conf. Acoustics, Speech, Signal Processing, pp. 2385-2388, 2001. 
[26] M. Al-Mahmoud, "Perfomance evaluation of code-spread OFDM," 46th Annual Allerton Conf, Allerton House, University 

of Illinois-Urbana Champaigne, Illinois, USA, pp. 274-278, Sept. 2008. 
[27] L. Hanzo, M. Miinster, B. -J. Choi, and T. Keller, OFDM and MC-CDMA for Broadband Multi-user Communications, 

WLANs and Broadcasting, John Wiley & Sons, 2003. 
[28] Z. Ye, G. J. Saulnier, D. Lee, and M. J. Medley, "FEC coding with a rate adaptive spread spectrum OFDM system," 

Communication Theory Mini-Conference, Vancouver, BC, Canada, pp. 67-71, June 1999. 
[29] M. G. Parker, "Close encounters with Boolean functions of three different kinds," 2nd International Castle Meeting on 

Coding Theory and Applications, Valladolid, September 2008. Also available at Lecture Notes in Computer Science (LNCS) 

vol. 5228, pp. 137-153, 2008. 

[30] F. J. MacWillams and N. J. A. Sloane, The Theory of Error Correcting Codes, Amsterdam, The Netherlands: North Holland, 
1986. 

[31] J. Seberry and M. Yamada, "Hadamard matrices, sequences, and block designs," Contemporary Design Theory: A Collection 
of Surveys, J. H. Dinitz and D. R. Stinson, Eds. John Wiley & Sons, Inc. 1992. 

[32] S. Litsyn, Peak Power Control in Multicarrier Communications, Cambridge University Press, 2007. 

[33] M. G. Parker, K. G. Paterson, and C. Tellambura, "Golay complementary sequences," Wiley Encyclopedia of Telecommu- 
nications, Edited by J. G. Proakis, Wiley Interscience, 2002. 

[34] M. G. Parker and C. Tellambura, "Golay-Davis-Jedwab complementary sequences and Rudin-Shapiro constructions," 
manuscript. 2001. Available in ' http://www.ii.uib.no/~matthew/ConstaBent2.pdf'. 

[35] S. B. Wicker, Error Control Systems for Digital Communication and Storage, Englewood Cliffs: Prentice-Hall, 1995. 

[36] I. S. Reed, "A class of multiple-error correcting codes and a decoding scheme," IEEE Trans. Inform. Theory, vol. 4, 
pp. 38-49, Sep. 1954. 

[37] J. H. Conway and N. J. A. Sloane, "Soft decoding techniques for codes and lattices, including the Golay codea and the 

Leech lattice," IEEE Trans. Inform. Theory, vol. IT-32, no. 1, pp. 41-50, Jan. 1986. 
[38] C. Tellambura, "Computation of the continuous-time PAR of an OFDM signal with BPSK subcarriers," IEEE Commun. 

Letters, vol. 5, no. 5, pp. 185-187, May 2001. 



